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Abstract 



We study the conversion of neutron matter into strange matter as a detonation 
wave. The detonation is assumed to originate from a central region in a spherically 
symmetric background of neutrons with a varying radial density distribution. We 
00 . present self-similar solutions for the propagation of detonation in static and collaps- 

^^ I ing backgrounds of neutron matter. The solutions are obtained in the framework 

f— ^ ■ of general relativistic hydrodynamics, and are relevant for the possible transition of 

^" . neutron into strange stars. Conditions for the formation of either bare or crusted 

strange stars are discussed. 



Oh; 

o 1 Introduction 

d '. Several authors have proposed that the ground state of matter is in the form of bulk quark 

matter made of the three flavors, up, down and strange. We adopt the term strange matter 
(SM) to refer to this bulk quark matter. Conversion of nuclear matter to SM is suppressed 



Vh ' at ordinary nuclear densities. At very large nuclear densities like those in neutron stars. 



where the Fermi energy is higher than the strange quark mass, the conversion however is 
possible. This has led to conjecture that strange stars that are predominantly made of 
SM may form by conversion of neutrons in dense neutron stars ^J |2]. The conversion is 
assumed to be triggered at the center of the neutron star (where the density is highest) 
either by presence of a seed SM or by spontaneous A -hyperons. 

Strange stars fall in a mass range similar to that of neutron stars, but they have 
smaller radiij^. There is mounting, albeit controversial, observational evidence for the 
existence of strange stars. Likely candidates are the the compact objects 3C58 and RX 
J1856.5-3754|ll|3]. 
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The dynamics of transition into SM has not yet been examined in details. According 
to El , the transition is reahzed by outward diffusion of strange quarks. The increased 
fraction of SM as a result of diffusion triggers further conversion into SM. In this case the 
conversion may be described as slow combustion at speeds of about lO^m/s. However, 
stability analysis in the presence of gravity has shown that this combustion process is 
unstable to small perturbations [?]. Once these perturbation grow they increase the 
effective area available for diffusion. Therefore the conversion may actually propagate in 
the form of detonation sustained by the enhanced "burning". In the simplest variant, the 
detonation wave (DW) is propagating into a neutron matter NM with constant density. 
Newtonian gravity was used in |H] to determine the parameters of the resultant strange 
star and to estimate the energy released in the process. 

More realistic models for the propagation of DWs should take into account a vary- 
ing density of NM. Further, next to black holes, neutron stars have the largest grav- 
itational fields in nature. Moreover, the characteristic time of gravitational processes 
t = (47rG^)~^/^ is comparable to the hydrodynamical time scale t = RNs/^d ~ 10~^sjn]. 
Therefore, it is prudent to consider detonation in the context of general relativistic hydro- 
dynamics. In this paper we study self-similar models of DW propagation in non-uniform 
NM, taking into account a self-consistent treatment of gravity according to general rela- 
tivity (GR). We consider detonation in static and collapsing backgrounds of NM. 

The paper is organized as follow. In §2 we write the equation of state for quark matter. 
In §3.1 we present the equations of special relativistic hydrodynamics in the absence of 
gravitation. Self-similar solutions in this case are presented in §3.2 for uniform and non- 
uniform backgrounds of NM. We turn to the general relativistic description in §4 where we 
derive the equations for general relativistic hydrodynamics and present the corresponding 
numerical self-similar solutions for static and collapsing non-uniform background. In §5 
we conclude with a discussion of the results and their implications. 

2 The Quark Matter Equation of State 

The equation of state (EOS) for hot (quark) matter has the form of simple relation j2], 

P= 3(^-^0), £0 = 45. (1) 

We use throughout e and p to denote the energy density and the pressure, respectively. 
This simple EOS may be obtained straightforward from the next EOS [H] 

P = ^^'T' + ^TV^ + i^/^' - ^' ^ = ^V + ^/i^ (2) 

e = g^^T^ + ^T^^ + ^^^' + B■, s = ^vr^T^ + 3T^^^ (3) 



T is the temperature, n is the baryon number density, /i is the chemical potential of 
quark matter, B is the MIT bag constant. The sound speed in such matter is defined 
as {dp/de)s = c^, where the derivation occurs at constant entropy. It is easy to see 
that the sound speed remains constant Cs = l/\/3. The EOS [T] was used by Witten for 
calculation of the structure of strange stars |2j. We shall adopt this simple approximation 
for modelling of detonation transition of neutron stars into strange one. 

3 Detonation in Special Relativistic Hydrodynamics 

First, we describe the propagation of detonation without taking into account the effects of 
gravity. This will be useful for comparison of our results with the existing literature and 
will also facilitate the interpretation of the general relativistic solutions to be presented 
in §11 The neglect of gravity is justified if the detonation time scale is smaller than the 
dynamical time scale. Gravitational effects become important after the end of detonation, 
but in this section we will only be concerned with describing the detonation phase. 

3.1 Basic equations 

Here we describe the relevant equations for detonation in special relativistic hydrody- 
namics (SRH). We assume that the pressure of NM is negligible compared to its energy 
density, 

Bn > Pn, (4) 

otherwise there will not be any significant gain in pressure in the transition from NM to 
SM. We will also assume throughout the energy density of the SM Eg behind the front is 
substantially larger than the "vacuum energy" e^ = 4:B which is close to nuclear energy 
density of usual matter, 

e, > 45. (5) 

Under the conditions(JH)-(jHl) we express these EOSs as 

p = ne , (6) 

where n = 1/3 and n = correspond to hot relativistic matter and cold neutron matter, 
respectively. Values other than n = can be used to describe approximately the EOS of 
the NM, in these cases n can be regarded representing the square of an "effective" sound 
speed for the NM. 

We will base our description of detonation in SRH on an analogy with the classical 
theory of detonation [H]. Therefore, instead of using covariant notation we express the 
relevant equations in terms of the radial component of the 3-velocity v{T,r) = dr/dr. 



where r is a radial coordinate and t = ct is the time measured in units of lengths. We 
use the notation w = e + p ior the "enthalpy" . 
The relevant equations of SRH are, 

1 / dp dp\ 1 f dv dv\ 
w \dr dr ) W^ \dT dr ) 
1 ( de de\ 1 ( dv dv\ 2v 



w ydr dr I W^ \dr dr I r 
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where W = Vl — 

The first of these equations is the Euler equation and the second is the continuity equation 

for energy density. 

These equations determine the evolution of the system in which a detonation front 
propagates outward from the center of the neutron star converting its into SM. In a frame 
of reference moving with the front, a requirement of the energy - momentum flow be 
conserved yields the following jump conditions at the front 

iPs + SsH/W! = {Pn + SM/Wh 

{Ps + SsMf/W^+Ps = {Pn + e^){vif/Wl + Pn. (8) 

In addition, there is baryon flux conservation n 

UsVi/W, = UnVi/Wn. (9) 

Here the indexes "n" , "s" designate the quantities concerning to NM and to SM, respec- 
tively. So, f^ is the velocity of "inflow" of neutron matter in front of DW and vl is the 
departure velocity of SM from the front. 

If the 3-velocity v^ in (jHJ is measured in the local frame of detonation front then, 
according to the relativistic equation of velocity transformation, in the laboratory local 
frame it takes the value 

v{D) = {D-v^)/{l-Dvf), (10) 

where D is the speed of DW in this frame. 

The classical analysis of detonation conditions leads to definition of a normal detona- 
tion condition when the speed of hot fluid corresponds to sound speed vl = Cs P]. In the 
case of NM at rest 

D = -v^ 

and substituting (jHl) in the jump conditions (jH)) yields 

v{ = cs: n = D{D-V?>/2)/{l-V?>/2D), e^ = D /{2/V?> - D)e^. (11) 



3.2 Self-similar solutions in SRH 



Let the detonation wave propagate in cold NM background (ra — > according to(^), and 
let the SM behind the front be hot and relativistic (see ©)• We will obtain self-similar 
solutions of ((Zj) for this situation. 

By analogy to classical detonation, we find the self-similar solutions assuming the 
velocity v to depend on r and r only through the variable 

^ = r/T. 

Analysis of the system ((7j) shows that self-similar solutions can be obtained in the case 
of power dependence of density of NM, 



iXQ/rf 



const. 



This distribution of NM can be created by its own gravity (Newtonian gravitation). But 
for a fast detonation wave (see (jllj) ). the inhomogeneity of NM distribution is only im- 
portant. (Similar situation is realized for a strong shock wave in an inhomogeneous 
atmosphere jTU]-) According to ()11|) . we have at the front position r = r^, 



3£:„ = 3 {ro/rdf e„. 



(12) 



and the front velocity remains constant. Therefore, the density of SM at the front varies 
with time according to the relation: 

Ss = 3(ro/r)''e„. 

This relation motivates the following form for v and e^, 

v = v{0, e = {ro/r)'e{0. (13) 

By expressing the partial derivatives d/dr and d/dr in terms of d/dC, the Eqs. ([7j) are 
transformed into the following ordinary differential equations. 



dv 




(^ — v) dv vk 



(14) 
(15) 



3.2.1 Detonation Wave in Homogeneous Background of NM 



We first present solutions for tlie propagation of DW in a liomogeneous distribution of 
NM (fc = 0). In the case the equations Eqs. (fT^ - (fTH|l become equations is 



dv 



2vW^ 



e(i-<) 
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d^ W^U-vUd^ 



and in the non-relativistic case {v,Cs << 1) the equation for v reduces to Zeldovich's 
classical equation for detonation 0, 



dv 
dC 



2v 



e ((v-O'-c^ 
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Figure 1: A phase plane of the solutions of (fT^ for propagation of detonation wave in a 
homogeneous (a) and an inhomogeneous background (b,c). In (a) there is a solution for 
velocity f (^) of gas that leaves detonation front propagating with velocity D and comes 
to node of solutions ^ = 1/a/3, v{^) = 0, where it may be matched to the solution for rest 
matter v = 0. For k = 1 (panel b) node shifts in the plan {C,,v). For k = 2 (panel c) 
there in no node but we have a unique solution leaving detonation front and converging 
tot; = 



Analysis of the first equation of p(j|l in the phase plane (^, v) reveals a variety of 
solutions which all converge to a node at ^ = c^ = 1/^/3 and v = with the derivative 
dv/d^ = 0. The solutions can be obtained independently for the ranges < ^ < Cg and 



the Cs < i < D. The solutions can be found by numerical integration of the equation (fTB|l 
starting from an arbitrary point in phase space. 

Examples of these solutions are shown in figure (0 a). As is clear from this figure the 
only physically accepted solution in the range < ^ < c^ is f(^) = (see [9 ) shown as 
the thick solid line. To the right of .^ = Cs the solution shown as the dashed dotted line 
corresponds to normal detonation. For this solution the velocity at the front S, = D is 
equal to the sound speed (see (fTTjl ) and is zero at ,^ = c^. The derivative of the velocity 
diverges as ^ approaches D, and tends to zero at ^ — > Cg. This detonation solution 
describes hot matter at rest for ^ < Cs and moves outward for Cg < ^ < D. The velocity 
of the front is \/3/2, according to (fTT|) . 

Our solution for the propagation of detonation in a uniform distribution of NM is 
fully consistent with previous calculations in the literature |HI. We also refer the reader 
to the literature treating the similar problem for relativistic detonation in "scalar field" 
stars [TT] . 

3.2.2 Detonation Wave in an Inhomogeneous Background of NM 

We search now for self-similar solutions of Eq. ()14|) for k > 0, i.e., for the propagation of 
detonation in NM with decreasing density. For k > 0, also, the equation has a node to 
which all solutions converge with the identical derivative. In this case, however, the node 
is at f > and not at f = as in the k = case. Examples of the solutions are shown in 
panels b and c of Fig.l. As in the homogeneous case, the physical solution corresponds 
to with f = at ^ = 0. The point ^ = f = is a critical one (saddle) and the solution 
originating from this point is the separatrix and is shown as solid line in the figure. Near 
^ = f = we find that 

v = ka^, (17) 

for ^ <^ 1. For k < 4/3 the node lies to the left of DW wave at ^ < D. Therefore, 
self-similar solutions corresponding to normal detonation (cf. [TT|) ) exist only for k < 4/3. 
For k > 4/3 the node shifts to the right of the DW. In this case, there are no self- 
similar solutions for normal detonation in a static neutron configuration. But we can find 
self-similar solutions for over-compressed detonation. For such type of detonation the 
velocity of the hot matter at the front is determined by the details of the actual physical 
situation but it is less than sound speed in SM jU] The position of DW front is defined 
from (jHl), f-^ = —Vd, 

Vd<Cs: D = —^-^ — -^; Ws = in + l)e„ 

{€s + Ps) Vd 

Then a laboratory velocity of outgoing SM from the front of wave ()10j) will be supersonic 
nearby the front of DW but transition through the sound speed in the self-similar solution 




occurs continuously and f ^ at ^ — > according to the formula (fT7|). 

The conditions (fTHj) are actually conditions of so-called "over-compressed" detonation. 

In our statement of the problem, at NM densities less than 10^^ g/cm^ we deal with a 
transition mechanism based on diffusion of strange quarks from SM to NM. Corresponding 
velocity of diffusion Vd depends on density of NM [HI, but in principle, it should be less 
than sound speed in SM. Therefore the velocity of SM leaving the front is equal to the 
velocity of the anomalous (turbulent) diffusion Vd-, which may be higher than the estimate 
given in |6j). 

4 Detonation in General Relativistic Hydrodynamics 

4.1 Basic Equations 

As we shall consider a spherically-symmetric problem, it is convenient to write the metric 
as 

ds^ = e^dr^ - e^dr^ - r^dil^; dil^ = dO^ + sin^ 6 dcp^ (19) 

Here the standard spherical system of coordinates r, 6, (p is used. Variables v, A are func- 
tions of r, T. 

To derive self-similar solution in form similar to SRH analysis (jTj) we will use an 
"auxiliary" time rj instead of the proper time of local observer r. The transformation 
from variables r, r to new variables 77, r is carried out with help of the ,Tacobian|12j: 

dit-r) dt (\-v\ 

The system of hydrodynamics equations in GR are then 

1 I dp dp\ 1 I dv dv\ 1 I dv d\\ , , 

w \dr drj J W^ \dri dr J 2 \dr drj J ' 

Iodine d\ne\ 1 f dv dv\ 2v 1 f dX diy\ 

\ dr] dr J W^ \dr drj J r 2 \d7] drj 

The Einstein equations for the metrics can be written down as: 

dX / 2\ ( (^^ ^^\ n ^^ K{p + e)rve^ 

drj ^ / yQj- Qj-J ' Qfj ]y2 

dX dX (e^-l A 

— + v— = -v\ + Kpre^\, (23) 

drj dr \ r J 



where G, c are Newtonian gravitation constant and speed of light, respectively and k = 
SvrG/c^. The equations PT|) - (PHj) can be presented as system of three equations by elimi- 
nating of function i^(r, rj), 



1 fOp dp\ I I dv dv' 
p + e \dr drj I W^ \ drj dr 



]_ dX dX 
2v \dr] dr 



' de de' 



1 



p + e \dri 



dr 



dv dv\ 2v V I dX dX 
yy2 yQj, Q^ j J. 2 \dr drj ^ 



dX 

drj 



k{p + e)rv e^ 



(24) 



In the limit of A = this system of equations reduces to the SRH equations (fT3j). 



4.2 Self-similar solutions in GRH 

In contrast to SRH, where self-similar solutions could be found for all k, in GRH the full 
system ()24|) permits only self-similar solutions only with NM distribution density ()12p 
with k = 2. We show in appendix El that this profile describes approximately a realistic 
density profile in a neutron star with a Harrison- Wheeler- Wakano equation of state |13j. 
In analogy with SRH we make the ansatz. 



^ = r/ri, Kp{ri,r) = 2n{^) /ri'^, v {ri,r) = v{^), A = A(0 



(25) 



Adopting an EOS of the form (jHl) and using v, U, X expressed in terms of the self-similar 
variable ^, the GRH equations are expressed as. 



dv 


' 2nvW^ 2nW^ 
[e(l-<) (n + l)(l- 


vO' 
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(26) 


d\nn ^ 


n + 1) 


\ h ( 2 
[{1-vO [n+1 


2v^ 


V {n + l)ne 


A " 
-/l 


//2 


(27) 


d^ -^ 


2(1 -O n 




dX 
dC 


2{n + l)ne^v 
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1 — f 4 


(28) 




nW^ 



In the case of hot SM {n = 1/3) , Eqs. (I2S1)-(I2H1) have a form. 



dv 



2vW'^ 



i{l-vi) 2(l-<)' 



?,W^ He^ I ( v-S 



l-< 



— V 






d\nU _ 4 ({^-v)dv V (v-^dX) dX _ Ue^v 

~d^~{l-vO[ Vr2 ^+2+ 2v 'd^y 'di~W^' ^ ' 

where 11 = 877. 

We first investigate the boundary conditions imposed on this system of equations (j^^ . 
We assume that at ^ = we have 

v{0) = 0, n(0), A(0) = const. (30) 

In the vicinity of ^ = we solution can be approximate, up to third order in ^, as 

v = a^, U = Uo-(3e, e-" = l + 7e^ (31) 

where 

a = 1/2, /3 = no(no/2-i), 7 = -no/4. 

4.2.1 Detonation in a static background of NM 

The static solution of the system ()28ll2(jj) for cold NM background 

^; = 0, i7 = Ci/e' e-^ = C2. (32) 

is also the analytical solution of the Oppenheimer-Volkoff equation (see Appendix)in the 
case n ^ 0, C2 = Ci{{n + l)/2n)^. 

This full solution of 123 for hot quark matter is found numerically and shown as the 
solid line in Fig|21 It corresponds to the separatrix by analogy to the same problem in 
SRH (shown as the dash-dot line in the figure) and is valid mathematically up to some 
^ < 1. However, the solution must be terminated at the DW position, D. In the case 
of normal detonation, the the position D on the axis ^ is found by the requirement that 
the speed of outgoing SM v-^ is equal to the sound speed. The previous analysis has 
shown, however, that for DW propagating in a background p„ ~ 1/r^ such point can 
not exist without taking into account gravitation. The gravitational influence results in 
backward motion of node towards the front of normal detonation in the phase plane but 
jump conditions of normal detonation exist only in hot neutron stars with n -^ 1/3. This 
dynamical process may convert hot NS into SS without any appreciable thermal energy 
release. In cold NS we must deal with over-compressed detonation (see section IH. 2. 2j) . A 
numerical calculation shows that v-^ can be a little less than Cg. 

For subsequent analysis of numerical models of transition we will be using analytical 
approximation of neutron star (seelHI). It satisfactorily describes distributions of density 
6, pressure 11 and the metrics A of neutron star from some initial dimensionless radius xq- 
Thus 

n = 8nA{n)/x^, A{n) = 2n/(n^ + 6n + 1), (33) 

-^11 = e-^ = 1 - 2A = const. (34) 
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Figure 2: The dynamical variables n(^), vgr{^) and e"'^*^^-* obtained from the self-similar 
solution of Eqs. ()29|1 . There is a small difference between the velocity in vgr and vsr 
obtained, respectively, from the general and special relativistic solutions. 



According to (fTHj) . the pressure of SM matter behind the DW front with account of 
will be defined by relation: 



4p(r/, r) = 4 n{Cj/7f = (n + l)T{D) A{n)/x^ 



(35) 



since the front of wave propagates with constant speed D and time r], as well as r, is 
measured in units of L* and t] = x/D. The initial moment of time is chosen so that the 
front of wave is found at some xq where the analytical approximation of distribution of 
NS density is true. (Note, that true distribution has a maximum in the center and a 
singular analytical solution approximates only the decreasing of this density with growth 
of X starting with finite Xq. It is supposed that the specified self-similar solution realizes 
from the moment of achievement of point xq by the DW front and further.) Thus for the 
representative of pressure at the front of wave we have 



U{D) = 2T{D){n + l)A{n)/D^' 



(36) 



And position of D is defined in the following way: at the point ^ = D values ^, v and 
v^ have to obey the law of energy-momentum flow conservation (jHJ and the relativistic 
equation of velocity transformation simultaneously: 



{2,{vff + l){n + l) ^n^ - 10^/3 + 1 



f\ 



1 + i{vi 



/(O (37) 



Thus, we are starting the numerical solution from point ^o = taking into consideration 
an analytical research of the solution in the vicinity of ^o fl30|) . At that, the setting of 



11 



Ho and Aq should be chosen so that after definition of the front position they would be 
corresponding to the mentioned values (jHUj) and (jMj) . 

4.2.2 Detonation in a collapsing background of neutron matter 

Recent observational data strongly suggest that "long-duration" Gamma Ray Bursts 
(GRBs) are associated with supernova explosions J14j. These GRBs follow the supernova 
events after a period varying from hours to months. This may point to the possibility that 
these GRBs appear as a result of the conversion of NM into SM during collapse of the SN 
remnant. This occurs when the mass of the NM becomes larger than the Chandrasekhar 
mass limit. The detonation in this case propagates on the background of the falling NM, 
so normal detonation wave forms. 

Negligible pressure forces in the collapsing NM allow us to neglect pressure in this 
statement of the problem. There are only two relevant parameters determining the self- 
similar solutions. These are the mass and the radius of the remnant strange star. The 
parameters of obtained self-similar solution are compare with those of the observational 
data. 




Figure 3: The dynamical variables obtained in the self-similar solutions in a collapsing 
background of NM with n = 0. 



In this statement of problem the jump conditions on the front (jHl) yield 

vi = 1/73, vl = V3/2, Ss = Ssn. (38) 

As in the stationary case, we are solving system (^^ for hot strange matter from point 
^0 to place of the front 



D 



1 — VVn 



(39) 
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The solutions of system ()26 p - (j^Hj) for NM and SM are presented in Fig.Q. At that we 
will claim the metric e"'^ = 1 — 2MG/ Rc^ to be defined by mass M of the strange star at 
the radius R{aX r = R pc^ = AB). Here we use parameters M = I.IMq and R = 9.7km 
(see data in 114| ). 



4.3 Numerical estimates 

As an example of detonation in static neutron star, let us consider the model of detonation 
transition of neutron star with mass of central part M ~ 0.59Mq and radius ~ 8.4 km 
measured in characteristic lengths L^, = 2.8 ■ 10^ cm that corresponds to x = 1.4. The 
analytical approximation of density distribution p„ measured in units of p* = 1.34 ■ 
lO^^g/cm^ is set by expression ^ with n = 1/12. At that A = 0.111, e'^ = 0.76. The 
results of numerical research for the over-compressed detonation are submitted in Fig. |21 

In the example under consideration the velocity of DW is equal ~0.85, that is close 
to the velocity D = ^3/2 = 0.851. 

It is interesting to estimate the value energy concerned with propagation of hot matter 
behind DW front, when the last one reaches the specified boundary radius of neutron 
configuration x = 1.4: 

-ft ix o 

E = An fr^ [t° - T^{v = 0)] dr = IGti f ^p r^dr (40) 



whence, using flS^ and (jH^ . we shall have 



The value of the dimensionless integral, obtained on the basis of the numerical solution, 
gives the value 0.057, corresponding to £" = 4.8 • lO^^erg = O.O27M0. This is close to the 
the energy produced in a gamma-ray burst p3|. 

The gravitational mass defect Am. According to [12j we have, 

R R 

Am = 47r /"t° (e^/^ - l) r^dr = An f p (— - l) (e^/^ _ ^j ^2^^ 


or, using (jH^j) and (|H^. 
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Thus,we have Am = 1.3 ■ lO^^g ^ o.OOTM©. 

The energy release associated with detonation can be roughly estimated as i? — Am. 
In our model it is only about 4% of the mass M. 

The detonation will stop when the NM density is close to AB/c^, leaving a crust 
made of the original matter of the neutron star. This crust can be expelled if the kinetic 
energy of the SM is large enough. The kinetic energy needed for that can be estimated 
by comparing the binding energy of the crust, 



Eb = 47r f {^pe^'^ - p^r'^dr 



R 

to the energy release, E — AM. We can see that in the case n = 1/12, Eh = O.OOSMq 
which is less than the energy release E — Am = O.OISMq. Therefore for this value n 
the crust might be expelled, leaving a bare strange star. Table 1 lists estimates for other 
values of n and for collapsing NM. 

Table 1: The parameters of detonation in static (f„ = 0) and collapsing (u„ = —0.13) 
background of neutron matter. Here pc is the central neutron density prior to detonation (in 
IQ-*^^ g/cm'^), M - mass of SM core, R is the radius at which detonation stops (in km), Vs -speed 
of SM behind the detonation front, E is the SM "kinetic energy", Am is E^ the binding energy 
of SM core and NM crust, correspondingly. All energies are in 10^^ erg. 



Pc 


n 


M/Mq 


R 


Vs 


E 


Am 


Eb 


9.40 


0.2 


0.59 


5.23 


-0.51 


35.18 


38.49 


20.94 


6.71 


0.16 


0.60 


5.42 


-0.52 


50.00 


21.35 


22.13 


5.37 


0.08 


0.573 


8.37 


-0.52 


48.44 


13.24 


28.27 


cohapse 





1.1 


9.74 


-i/Vs 


87.72 


58.64 


- 



5 Discussion 

We have studied the conversion of neutron matter into strange matter as a detonation 
wave propagating outward of a central region. The dynamics of neutron-quark transition 
has been studied earlier for a uniform static background of neutron matter in the context 
of special relativity |n|. Here we employ special and general relativistic hydrodynamics to 
find self-similar solutions for the propagation of detonation in a non-uniform background 
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of neutron matter. Both static and collapsing backgrounds are considered. Solutions 
for detonation in a static background may correspond to the conversion of neutrons in 
an otherwise stable neutron star. In a collapsing background the solutions may describe 
conversion in the neutron core of SN Type II ^1 ^^. In this latter case detonation is 
assumed to be triggered at the onset of the collapse. 

A description of the physical processes leading to detonation can be found in pa- 
pers jni U\ ■ Briefly, the initial stages of conversion of neutrons to strange quark matter 
proceeds by diffusion of SM into the neutron matter. The diffusing SM then stimulates 
further conversion and the whole process may be described as combustion [^. However, 
the process turned out to be hydrodynamically unstable when gravity is taken into ac- 
count J7j. Hydrodynamical instabilities increase the surface area available for diffusion 
allowing for the possibility of propagation in the form of detonation. In the other words, 
these processes leads to anomalous diffusion of strange quarks to the detonation front 
with speed that is equal to speed of leaving hot matter. 

In special relativity the self-similar solutions describe detonation in a static background 
of NM with distribution p ~ r^^. It was shown the condition of a normal detonation exist 
for k < 4/3, while denotation must be over-compressed for k > 4/3. 

Self-similar solutions of the problem in general relativity are permitted only for k = 2. 
The reason for this restriction is the appearance of new fundamental constant G in GR. 
This additional physical scale prevents the development of a self-similar behavior of the 
type considered here. In special relativity the detonation for k = 2 is over-compressed 
independently of the equation of state of cold matter. However, in general relativity we 
obtain normal detonation for n close to 1/3 (P = ne). Below a limiting value of n only 
over-compressed detonation is possible. We have not been able to find the limiting value 
of n analytically but our numerical calculations show that it should be n ~ 0.22. 

Normal detonation is obtained for conversion in collapsing cold NM. In this case the 
pressure of NM is completely neglected, i.e., n = in (jH)). In such statement of problem 
one can easily calculate all relevant quantities using only two parameters characterizing 
the system, e.g., the mass and radius of the strange star. 

The propagation of detonation in neutron stars must halt at the radius where the 
neutron density is close to the quark vacuum energy, p = AB/c. When the DW reaches 
this radius the SM would be surrounded by a neutron crust. The description in this pa- 
per breaks down at this stage of the evolution. However, we still can use the self-similar 
solutions to determine the energetics of the system and to predict some general features 
its subsequent evolution. In all of our self-similar solutions, except those describing det- 
onation in a uniform background, the velocity of SM is moving outward. Although the 
"kinetic energy" of the SM is insufficient to disassemble the whole configuration (see Ta- 
ble 1), it might be enough to expel the crust surrounding the quark core. According to 
Table 1, if n is close to 1/3, the crust is maintained. In this case we may end up with 
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a crusted strange star. For lower n the energy is sufficient to remove the crust, leaving 
a a bare strange star. In our estimates presented in Table 1 we have neglected a few 
factors. Non-vanishing mass of strange quark and constant of strong interaction taken 
into account will lead to reducing of energy release. We also have neglected any further 
conversion (e.g., by slow combustion in the crust) after the detonation is halted. 
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6 Appendix A: Density profile of Neutron stars with 
a HWW equation of state 

In this section we solve the Oppenheimer-Volkoff equation to compute the density profile 
in neutron stars with the Harrison- Wheeler- Wakano (HWW) equation of state 

1 / (5^/3 \ 
^ " 3 1^(1 + 55/9)6/5 j ' (42) 

where H = p/p^c^, 6 = e/p^c^, p, = 1.34 ■ lO^^g/cm^. The HWW EOS is valid for 
central densities > 4.63 ■ lO^^g/cm^ ^2|- The sound speed c^ is given by 



(43) 



2lV _dp_dll_ 1 6^/^(5 + 36^/^) 
~) ~Te~~di^~9 (1 + (55/9)11/5 • 

Therefore for small 6 this EOS describes cold NM. 
From equations (jSHl), after integration, we have 

r 

e"^ = 1 - {2G/c^)m{r)/r, m{r) = Air/c^ e{r)r'^dr , (44) 



where e{r) is the energy density of NM and m{r) is the mass of sphere within a radius r. 
The Oppenheimer-Volkoff for the equilibrium of NS is fH] 

dp {e + p) G[m{r) + 4:'Kr^p{r)/c^] , , 

dr rc'^[r — 2(G'/c^)m(r)] 
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To present Eqs. ()44|) - ()45|) in dimensionless form, we use the radial coordinate x defined as 
X = r/L^ where the characteristic length L^, = cj ^^.i^p^G = 2.8 ■ 10^ cm. In dimensionless 
units the mass within a radius r is 



/i(x) 



5{x)x^dx. 



(46) 



In order to get the physical mass we multiply fi by m^ = AnL'^p^ = 3.8 ■ 10 g = 1.91Mq. 

Let us denote p*Ll = 2/kc^, and if arguments {r,ri) are measured in units of L^,, 
dimensionless pressure 11 correlates with representative of pressure -/I(0 as 11 = n{^)/iif. 

Oppenheimer-Volkoff' s equation ()45|) in dimensionless variables can be presented as 



dn (5 + n)(nx3 + /i(x)) 



dx 



x{x — 2yu(x)) 



(47) 



This equation together with Eq. (j46|) and EOS (j42|) completely defines the structure of 
the neutron star. Given the density at the center and assuming a vanishing pressure 
at the star's boundary, a unique solution to these equations can be found by numerical 
integration. We denote the density at the center (x = 0) by 5o. The numerical solutions 
are shown as the dashed lines in the three panels in FigEJ Also, the solid lines in these 
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Figure 4: The profiles of the density distribution 5, pressure H, and mass fi as function 
X. The dashed hues are obtained for the HWW EOS with central density of 6 = 0.4. The 
solid lines correspond to the analytical approximation (0S)) with n = 1/12. 



panels show solutions obtained by assuming 

5 = A/x\ U = n6 



(45 
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Thus we have: 

n = An/x^] /i(x) = Ax. (49) 

It is easy to see that this special solution satisfies to Eq. ()47p where the coefficient A is 
defined by expression: 

A = 2n/{ri^ + Qn+l). (50) 

At that the components of the metric {—g^^) and f^oo = e^ are given by the formula 

e-^ = l-2A, u = B{n)\Yi{x/x.), B{n) = '^^^^^, (51) 

where x < x^, (here x* is an integration constant). It is interesting to note that the {—g^^) 
is constant, while g^Q depends on x. 

This solution has a conic singularity at x = 0. But the metric at this point has no 
singularity. We approximate density profiles of neutron stars by shifting this solution at 
X ^ X + xq. Comparison of this simple analytical solution with a numerical calculation 
is presented in Fig. HJ 

Therefore approximating the NM distribution with a power law oik = 2 as satisfactory 
away from the center. 
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